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Abstract 

In the present work, the temperature dependence of the scalar mesons parameters 
is investigated in the framework of thermal QCD sum rules. We calculate a -pole and 
the non-resonant two-pion continuum contributions to the spectral density. Taking into 
account additional operators appearing at finite temperature, the thermal QCD sum rules 
are derived. The temperature dependence of the shifts in the mass and leptonic decay 
constant of scalar cr(600) meson is calculated. 

1 Introduction 

The QCD sum rules method p], proposed about three decades ago, is one of the powerful 
methods for investigating the properties of hadrons. This method has been extensively used as 
an efficient tool to study the masses, decay form factors and so on [2]. In recent years there 
has been increasing interest in the modification of hadronic properties at finite temperature in 
order to understand the results of the heavy ion collision experiments. 

The QCD sum rules method is extended to the finite temperature in [3] and finite tempera- 
ture sum rules have several new features. One of them is the interaction of the current with the 
particles of the medium. This effect requires modifying hadron spectral function. The other 
novel feature is the breakdown of Lorentz invariance by the choice of reference frame [l]-[6]. 
Due to the residual 0(3) symmetry more operators with the same dimension appear in the op- 
erator product expansion (OPE) at finite temperature compared to those at zero temperature. 
Taking into account both complications, the investigation of OPE for thermal correlator of the 
two vector currents, and the thermal QCD sum rules for vector mesons have been realized in 
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[7] and [8], respectively. Also, nuclear medium modifications of meson parameters are widely 
discussed in the literature [9]-[T0]. 

In the present work, we investigate the properties of the scalar a meson in the framework 
of thermal QCD sum rules. The basic idea of thermal QCD sum rules is to get information 
about temperature dependence of hadron parameters, studying the same correlator, both at 
high temperature where the quark-gluon plasma is expected and at low temperature, where the 
hadronic phase is dominated. Note that the nature of light scalar mesons is still an open problem 
and is the subject of intensive and continuous theoretical [TT] and experimental investigations 
[12J. Can we get any new information about the nature of the scalar mesons from the thermal 
QCD analysis? Present work is addressed to the investigation of this problem. 

The paper is organized as follows. In section 2 we derive the thermal QCD sum rules for 
scalar er(600) meson. In section 3 we present our numerical calculations. This section also 
contains discussion and our conclusion. 

2 Thermal QCD sum rules for scalar sigma mesons 

In this section we construct the thermal sum rules for scalar er(600) meson. For this purpose 
we consider the thermal average of correlation function 

T(q) = i J d 4 xe tq - x (T(J(x)J(0))), (1) 

where J(x) = -^(uu + dd) is the interpolating current with the a meson quantum numbers. 
The thermal average of any operator is determined by following expression 

(O) = Tre- pH 0/Tre~ pH , (2) 

where H is the QCD Hamiltonian, and (3 = 1/T stands for the inverse of the temperature T 
and traces are carried out over any complete set of states. 

The fundamental assumption of Wilson expansion is that the product of operators at differ- 
ent points can be expanded as the sum of local operators with momentum dependent coefficients 
in the form: 

T( q )=Y,C n (q 2 )(O n ), (3) 

n 

where C n (q 2 ) are called Wilson coefficients, and O n are a set of local operators. In this expan- 
sion, the operators are ordered according to their dimension d . The lowest dimension operator 
with d = is the unit operator associated with the perturbative contribution. In the vacuum 
sum rules operators with dimensions d = 3 and d — 4 composed of quark and gluon fields 
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are the quark condensate (ipip) and the gluon condensate {G a G atMU ), respectively. At finite 
temperature Lorentz invariance is broken by the choice of a preferred frame of reference, and 
therefore new operators appear in the Wilson expansion. In order to restore Lorentz invariance 
in thermal field theory, four-vector velocity of the medium is introduced. Using four-vector 
velocity and quark/gluon fields, we can construct a new set of low dimension operators (uQ^u) 
and (uQ 9 u) with dimension d = 4 , where Q* and Q 9 are fermionic and gluonic parts of 
energy momentum tensor O^, respectively. So, we can write thermal correlation function in 
terms of operators up to dimension four: 

T{q) = CJ + C 2 (^) + C 3 (G; u G a n + C A (uQ f u) + C 5 (uG 9 u). (4) 

The Wilson coefficients in Eq.fllJ) are calculated in [13], and renormalization group improved 
expression of OPE leads to the following result 



, 4 Aju-Q) 2 1 
16 + 371/ V Q 4 Q 2 



1 6 

(uQu) + A (Q 2 ) ( y (ue f u) - («9%) 



/ 

2M~<V 6 



(5) 



where Q is the Euclidean momentum, \(Q 2 ) = [a s (fi 2 ) / a s (Q 2 )] and Mi , = n/Qf + 0^,. 
At one loop level the constants S and b are given by 

2/16 \ 2 
5 = - ( y + n/ J and 6 = 11 - -n/ , (6) 

where n/ is quark flavours number. The spectral representation for the correlation function in 
go at fixed |q| can be written as [8] 



/ \ f°° #N (q'o, |q|) 

T [Qn, q| ) = / ^Vo — 5 + subtraction terms 



(7) 



where 



AT(g ,|q|) = -JmT (g , |q|)tanh(/%,/2) and Qo = ~9o- (8) 

7T 

Note that, the subtraction terms are removed by the Borel transformation. For this reason 
in further discussion we omit these terms. Equating spectral representation and OPE, and 
performing Borel transformations with respect to Ql from both sides Eq.(JZJ), we obtain the 
QCD sum rules 



POO 

/ dg V^ 2 iV( go ,| q |) = e -iqi : 
Jo 
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where M is Borel parameter, (Oi) and (O2) are the non-perturbative contributions of higher 
dimensional operators, 



g 2 



(d) = 3m(#) + ^(G« u G a n , (10) 



<o 2 > = 12 



32vr 

-16 



(m9m) + A(M 2 ) ( — (ite'it) - ( M 9 9 «)j 



lQ + 3n f 

Now we consider the phenomenological part of the correlation function. We shall work below 
the critical temperature, where the physical spectrum is saturated by hadrons . In this case, 
similar to the vacuum QCD sum rules, the dominant contribution to the spectral density comes 
from u mesons. We also calculate the contribution of the non-resonant two-pion continuum. 

Let us calculate a-pole contribution to the correlator. The leptonic decay constant A CT of 
the cr-meson is given by, (0| J(0)|er) = m a X a , where m a is cr-meson mass. In thermal field 
theory, the parameters m a and A CT must be replaced by their temperature dependent values. 
The vacuum value of the leptonic decay constant is obtained from two point QCD sum rules 
and \ a = 200MeV [H]. The absorptive part of the correlation function is calculated by using 
the following field-current identity 

J(x) = m a X a a(x) , (12) 

and cr-meson contribution to thermal correlator can be written as 

T{q)=im 2 XDl x {q) . (13) 

Here D xl (q) = J d A xe iq ' x (T(a(x)a(O))) is the time ordered product of two cr-meson fields (11- 
component of the finite temperature scalar field propagator with mass m a in the real time 
formalism) and has the following form [T5]-|16j 



DIM = q 2_nl + i£ + 2 ™^ 5 ^ - O > ( 14 ) 



where n(u g ) is the Bose distribution function, n(u> q ) = [exp((3uj q ) — l]^ 1 and u q = ^q 2 + ra£. 
The imaginary part of correlation function can be simply evaluated using the formula = 
tt5(x) + iP which leads to 

ImT{q) = nml\l(2n(u q ) + 1)5 (q 2 - m 2 a ) . (15) 

With the help of 5-function we obtain the following result for cr-pole contribution to the spectral 
function 

N(q)=ml\l5(q 2 -ml). (16) 



In order to calculate the dependence of m a and \ a on temperature, we consider appropriate 
loop diagrams. Let us calculate the 7T7r-contribution to the amplitudes, which describes the 
interaction of the current with the particles in the medium. This contribution to the correlation 
function can be written as 

T(q) = igl J -0^D- u (k)DUk - q) , (17) 

where D^ik) is the 11-component of the finite temperature propagator for pions and g a = 
2,0 GeV [H]-[T8]. The integration over k in Eq. (fl7|) can be evaluated using the residue the- 
orem. After integration and some simplifications for the imaginary part of the correlation 
function we obtain 

f dk 1 

ImT(q) = irgl J —-^ — - ((1 + n x ) (1 + n 2 ) + rnn 2 ) (5 (q - u x - uj 2 ) + 5 (q + ut + u 2 )) 
+ ((1 + ni) n 2 + (1 + n 2 )ni) (5(q - lo x + u 2 ) + 5 (q + u x - io 2 )) , (18) 

where, 

ni = nM , n 2 = n(u 2 ) , lq x = \Jk 2 + ml , tu 2 = y/(k- q) 2 + ml . (19) 

At values q = LO\ + u 2 and q = ui\ — uj 2 the terms involving the density distributions can be 
written as 

[(1 + m) (1 + n 2 ) + nin 2 ] tanh (^j^j = (m + n 2 + 1) , (20) 

[(1 + m) n 2 + (l + n 2 ) m] tanh (^-j = (n 2 - m) , (21) 

respectively. As can be seen, delta function 5(q — uj\ —oj 2 ) in Eq. (fl8|) gives the first branch cut, 
q 2 > Ami, which coincides with zero temperature cut that describes the standard threshold 
for particle decays. On the other hand, delta function <5(g — + oj 2 ) in Eq.(fl8l) shows that 
an additional branch cut arises at finite temperature, q 2 < , which corresponds to particle 
absorption from the medium. Therefore, delta functions 5(qo — u)\ — uj 2 ) and 5(qo — <jJ\ + uj 2 ) 
in Eq.( fl8l) contribute in regions q 2 > 4ml an( i q 2 — 0> respectively Taking into account both 
contributions, the spectral function can be written as 

N(q) = gl J ^f^ 9 [("i + n 2 + 1) g (g 2 - iml) + (n 2 - m) (-q 2 )' 

x5 (q 2 - 2q uj l + 2|k||q| cos^) . (22) 

The integration over angle 9 in Eq. (j22l can be evaluated using the constraint | cos#qkJ < 1 , 
which leads to following inequality 

l„2 , , I 

< 1 • (23) 
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The solution of this inequality at values q 2 > Am 2 give us the integration range of U\ as 
c<j_ < u) < lo + , where 

(24) 



w± = ^(q ± |q|u) 



v(q 2 



1 — Ami I Q 2 



(25) 



At q 2 < the region of variation of U\ must be oj + < u)\ < oo. Finally, the thermal spectral 
function can be written as 



N(q) 



gl /•<"+ dux 



r w? ( - + + :> 6 & - 4m -) + If c w? (n " - ni) ° ^ ■ <26) 



2|q| 



Changing the variable uj\ to x given by uj\ = | (go + |q|#) > we finally get the two pion contri- 
bution to the spectral function as 



N(q) ^ gf-^l + Nl (q) 



9lv(q 2 



+ 



8vr 2 J-, 



dxn(-(q + \q_\x] 



q 2 > Am 2 (27) 



N(q) = N 2 (q) 



16tt 2 



dx 



n Q (|qk - 9o)) - n Q (|q|« + ?o)) 



r < o 



(28) 



The QCD sum rules are obtained by equating theoretical and phenomenological parts of 
correlation function. Taking into account of expressions Eq. (l27|) and Eq.(28) in Eq.Q we get 

ml(T)\ 2 a (T)e~ m2 ^/ M2 + ^-M 2 ^ f°° ^V^^g 2 ) 

on J4m2+|q|2 

iqi 2 



|q| 2 /Af 2 



+ e 
3M 4 



4m2 + q|5 



<Oi> 



dqle-^N, (g , |q|) + jf ^Vo/^ | q |) 



(29) 



As the temperature approaches to zero, the two terms in bracket go to zero and the thermal 
average of the operators on the right become the expectation values, recovering the vacuum 
sum rules [H]. In the limit |q| — > 0, the sum rule (29) simplifies considerably. Finally we obtain 



that 



where 



ml(T)\l(T)exp(-ml(T)/M 2 ) + J (M 2 ) + h(M 2 
I (M 2 ) 



9l < : " 



h{M 2 
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07T Z J 4m 2 

2 coo 



3M 4 
" l^ 2 

ds v(s)exp(—s/M 2 ) , 
ds v(s) n(yfs/2)exp(-s/M 2 ) , 



+ (0) 



An 2 J Ami 

(O) = ((h) + (o 2 ) 



(30) 

(31) 

(32) 
(33) 
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3 Numerical analysis of the shifts in mass and leptonic 
decay constant 



In this section we present our results for the temperature dependence of the shifts in a meson 
mass and leptonic decay constant. By derivativing with respect to 1/M 2 from both sides of the 
sum rules (30), and making some transformations we obtain 

m 4 a Xlexp(-ml/M 2 ) - J X (M 2 ) + rj{Os} 



<{T) 



ml\lexp(-ml/M 2 ) - h(M 2 ) + (O) 



A=(T) = \] 



m 2 X + {(0) - h(M 2 )) 


exp(m 2 /M 2 ) 


™IK + O - i) 




J X (M 2 ) -r](0 3 ) +m 


1 ((0) - h(M 2 )) 


expim 2 / M 2 ) 



(34) 
(35) 



where the bar on the operators means subtractions of their vacuum expectation values and 

SM 2 



r](M 2 



Mn(M 2 /A 2 



Ji(M 2 



(Os> 



All 2 J Ami 
12 



ds s v{s) n (y/s/2) exp(-s/M 2 ) 
16 



-\(M A 



(uQ f u) - (uQ 9 u) 



(36) 

(37) 
(38) 



16 + 3n f v ' V 3 

For the numerical analysis, let us list thermal average of operators contributing to the QCD 
sum rules. The temperature dependence of quark condensate is known from chiral perturbation 
theory p2]-[20] 



<^> = (o\W\o) 



1 - 



nj-1 T 2 



+ 0(T 4 ) 



(39) 



n f 12F 2 

where nj is number of quark flavors and F = 0.088 GeV. The low temperature expansion of 
the gluon condensate has been studied in article 



g_ 

At: 2 



(40) 



/ / 

where the trace of the total energy momentum tensor is given by (O^) = (@) — 3p, and for 
two massless quarks in the low temperature chiral perturbation limit the trace has following 
form [19 

— 0{T«>) . (11) 



270 T ; 

Here (0) is the total energy density and p is the pressure, whose expressions are known in 
the low temperature region [20]. The pion decay constant has the value of F n = 0.093GeV^ 
and the logarithmic scale factor is A p = 0.275GeV. We also use the fact that, the quark and 



gluon energy densities at finite temperature can be expressed as rif{Qf) = (Q g ) = |(0), which 
agrees both with the naive counting of the degrees of freedom and empirical studies of the pion 
structure functions [1], [H]. 
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Figure 1: Shift in the sigma meson mass as a function of temperature for M 2 = 1 GeV 2 and 
M 2 = 1.2 GeV 2 . 

For the numerical evolution of the above sum rule, we use the values miipip) = —0.82 x 
10~ 4 GeV A , A = 0.230 GeV and m a = 0.6 GeV. We study the dependence of a meson mass 
and leptonic decay constant on M 2 , when M 2 changes between 0.9 GeV 2 and 1.4 GeV 2 . This 
region of M 2 is obtained from the mass sum rule analysis of the a meson |14j . 

The shifts in a meson mass and leptonic decay constant as a function of temperature 
for different values of M 2 is shown in Fig.l and Fig.2, respectively. As seen, the results for 
Arricr are stable for temperatures up to 120 MeV. At high temperatures the results for Am a 
becomes unstable and the contributions of higher dimensional operators become important 
here, whose inclusions might restore the stability in M 2 to higher temperatures. The results 
for AA CT are stable and leptonic decay constant decreases with increasing temperature and 
vanishes approximately at temperature T = 160 MeV . This situation may be interpreted as 
a signal for deconfinement and agrees with heavy-light mesons investigations [23]. Numerical 
analysis shows that the temperature dependence of AA CT is the same, when M 2 changes between 
0.9 GeV 2 and 1.4 GeV 2 . 

Obtained results can be used for interpretation heavy ion collision experiments. It is also 
essential to compare these results with other model calculations. We believe these studies to 



8 




T (MeV) 

Figure 2: Shift in the leptonic decay constant as a function of temperature for M 2 = 1.2 GeV 2 . 

be of great importance for understanding phenomenological and theoretical aspects of thermal 
QCD. 
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